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Abstract 

The b — > ssd and b — > dds decays are highly suppressed in the SM, and are thus good probes of 
new physics (NP) effects. We discuss in detail the structure of the relevant SM effective Hamiltonian 
pointing out the presence of nonlocal contributions which can be about \- 4 (m 2 c /m 2 ) ~ 30% of the 
local operators (A = 0.21 is the Cabibbo angle). The matrix elements of the local operators are 
computed with little hadronic uncertainty by relating them through flavor SU(3) to the observed 
AS = decays. We identify a general NP mechanism which can lead to the branching fractions of 
the b — ► ssd modes at or just below the present experimental bounds, while satisfying the bounds 
from K — K and Bi s \ — Br s \ mixing. It involves the exchange of a NP field carrying a conserved 
charge, broken only by its flavor couplings. The size of branching fractions within MFV, NMFV 
and general flavor violating NP are also predicted. We show that in the future energy scales higher 
than 10 3 TeV could be probed without hadronic uncertainties even for b — > s and b — > d transitions, 
if enough statistics becomes available. 
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I. INTRODUCTION 



The decays b — > ssd and b — > dds are highly suppressed in the SM: they are both loop 

ements V t . and/or V t d). As such they 



3,|a|. The types of NP 



and CKM suppressed (by six powers of small CKM e 
can be used for searches of New Physics (NP) signals l|, |2|, |3j, |4], |5|, |6 
that would generate b — *> ssd and b — > dds transitions will commonly also give contributions 
to K — K, B — B and B s — B s mixing. Since no clear deviations from the SM predictions are 
seen in the meson mixing, is it possible to have deviations in b — > ssd and b — > dds transition 
observable at Belle II and at LHCb? A related question is: with improved statistics, can the 
experiments using b — > ssd and b — > dds decays push the bounds on flavor violation scale 
beyond what can be achieved from the mixing observables? 

We address the second question first. For simplicity let us assume that NP contributions 
can be matched onto the SM operator basis, so that H AS = C{ d {di, r y^si){di, r y l j,s£) , H AB = 
Q s (s L ^b L )(s L j^b L ) and H^ ssd = C^ ssi (s L ^b L )(s L -?^d L ). Using \C\\ = 1/(A*) 2 one finds 



B(B° - R»K») = 0.3 x 10" 6 (^^) 4 , (2) 



K — K mixing : A sa > 1.0 • 10 d TeV, 
B d - B d mixing : A 6d > 210 TeV, (1) 
B s - B s mixing : A bs > 30 TeV, 

with Im(C± d ) additionally constrained from ek- The above bounds should be compared with 
the following prediction for the b —>■ ssd transition in the presence of NP with scale A b ^ ssd 
(see section \V\ for derivation) 

10 TeV\ 1 

fab— >ssd 

while the SM prediction for this branching ratio is of O(10~ 15 ). Let us take as an estimate 
^b^ssd ^ y/ A 6s A sd , a relation that holds in a wide set of NP models including the Minimal 
Flavor Violation (MFV) and Next-to-Minimal Flavor Violation (NMFV) frameworks. With 
enough statistics the bound on A bs can then be pushed up to 10 3 TeV and higher without 
running into SM background. The b — > ssd decay modes could thus be used to constrain the 
NP flavor structure for b — * s transitions as precisely as it is possible for s — > d transitions 
from kaon physics. However, the statistics needed is very large. For instance, even to 
probe this type of flavor violating NP beyond the mixing bounds, the LHCb and Belle II 
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luminosities will not be enough. In this scenario the K — K and Bd — Bd mixing bounds 
translate to B(b — ► dds) < 1CT 13 and the bounds from K — K and B s — B s mixing translate 
to B{b -> ssd) < KT 11 . 

Does this mean that any NP discoveries using 6 — > dds and 6 — > ssd transitions are 
excluded at Belle II and LHCb? Certainly not. It is possible to have significant effects in 
b — > dds and b — > ssd while obeying the bounds from the meson mixing, if (i) the exchanged 
particle (or a set of particles) X carries an approximately conserved global charge and, if 
(ii) additionaly there is some hierarchy in the couplings (or alternatively some cancellations 
in K — K mixing). Consider the NP Lagrangian of a generic form 

^fl av or = g b ^ s {sTb)X + g s ^ b {bTs)X + g d ^ s (sTd)X + g s _ d (dTs)X + h.c, (3) 

and assume that X carries a conserved quantum number broken only by the above terms. We 
also assume for simplicity that the field X couples to a fixed Dirac structure T. Integrating 
out the field X produces flavor- changing operators 

£ cff = T^[^ s #^d(srd)(sf d) + g b ^ s g* s ^ b (sTb)(sTb) 

M x (4) 

+ 9b^ s g* s ^ d (sTb)(sTd) + g d ^ s g* s ^ b (sTb)(sTd)] , 

with the terms in the first line contributing to K — K mixing and B s — B s mixing, and in 
the second line to b — > ssd decays (we also introduced F = 7 r"''7 ). It is now possible to set 
contributions to meson mixing to zero, while keeping b — > ssd unbounded. This happens for 
instance, if 

9b-*a < 9s^b, g s ^d < 9dr-*s, or 9b-* s > 9s-+b, 9s^d > 9d^s- (5) 

In this way all the present experimental bounds can be satisfied, while branching ratios for 
b — > ssd and b — > dds induced decays are O(10 -6 ) (see section M for details). 

The important ingredient in the above argument was that X carried a conserved quantum 
number, so that there were no terms in £ e ff of the form 

gb U s (sTb)(sTb) + g 2 d ^ s (sTd)(sTd) + g* s l b (sfb)(sfb) + gf^sT d)(sf d) . . . , (6) 

These would be generated for X = X\ which is impossible, if X carries a conserved charge. 
If terms ([6]) are present, then B s — B s mixing forces both g b -> s and g s ^b to be small, and 
the hierarchy in (jSJ) is not possible (similarly K — K mixing bounds gd-* s and g s ^d to both 
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be small). An explicit example of a NP scenario where only terms of the form (J3J) are gen- 
erated is R— parity violating MSSM [4j. The i?-parity violating term in the superpotential, 
^ = Kjk-LiQjdk, leads to ViqLjdkR flavor violating coupling. Sneutrino exchange gener- 
ates operators of the form (jl]), while operators of the form (jH]) are not generated, since the 
sneutrino carries lepton charge broken only by .R-parity violating terms. 

A hierarchy of couplings in dSJ) is also present in (N)MFV models, if left-right terms give 



dominant contributions 10j. Both terms in and ([6]) are generated, on the other hand, 
for FCNCs induced by Z' exchange, since Z' does not carry any conserved charge. 

In this paper we will not confine ourselves to a particular model but keep the analysis 
completely general using effective field theory. We will improve on the existing SM predic- 
tions, and also give predictions for general NP contributions. The most general local NP 



hamiltonian for b — > ssd transition is 



2] 
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2~(J2 c jQj + Yl ^Qj) > ( 7 ) 



3=1 3=1 



A? 

where Cj are dimensionless Wilson coefficients, Anp the NP scale, and the operators are 

Qi ={sLl^b L )(s L Yd L ), 

Q2=(sRb L )(sRd L ), Qs = (s a R b P L )(4dl), (») 
Qa =(s R b L )(s L d R ), Q 5 = {s a R bi){sld a R ). 

The Qj operators are obtained from Qj by L <->• R exchange. In SM only Qi is present. The 
b — > dds effective Hamiltonian is obtained by exchanging s d in the above equations, while 
the K — K and B s — B s mixing Hamiltonians follow from b — > d and d — > s replacements. 

The predictions following from the NP hamiltonian require calculating the QCD 
matrix elements of the four-quark operators. In this paper we show that the Q\ matrix 
elements can be related by SU(3) flavor symmetry to linear combinations of observable 
AS = decay amplitudes. This gives clean predictions for the branching fractions of the 
exclusive b —>■ ssd and b — > dds modes in the SM and the NP models where Q\ dominate. 
This happens in a large class of NP models, including the two-Higgs doublet model with 
small tan/3, and the MSSM with conserved R parity [4j. The effects of the operators with 
non-standard chirality can be estimated using factorization. 

The outline of the paper is as follows. In Section [III we review the structure of the 
effective Hamiltonian mediating the b — > ssd, dds decays in the Standard Model. We point 
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FIG. 1: Matching the box diagrams contributing to b — > ssd decays onto an effective theory with 
in>W > A 4 — m b- The top quark box diagram (above) is matched onto a local four-quark operator, 
while the box diagrams with u, c internal quarks (below) are matched onto local and nonlocal 
operators. The mixed top-charm and top-up box diagrams are power suppressed by w? u Jm^ and 
do not contribute at leading order, see Appendix A. 



out that in addition to the local operators, the effective Hamiltonian contains also nonlocal 
operators which have not been included in the previous literature. In Section HTT1 we derive the 
flavor SU(3) relations for the matrix elements of the Q\ operator. The resulting numerical 
predictions for b — ► ssd, dds decays in the SM are given in Section IIVI NP predictions 
in the case of Qi operator dominance are discussed in Section |V] while in Section IVII the 
modifications needed for a general chiral structure are given. Three appendices contain 
further technical details. 



II. SM EFFECTIVE HAMILTONIAN FOR b -» ssd AND b -> dds DECAYS 



In the SM the b — > ssd, dds decays are mediated by the box diagram with internal u, c, t 
quarks, Fig. [TJ For notational simplicity let us focus on the case of b — > ssd, while the results 



for b — ► dds can be obtained through a replacement s <-> d. The effective weak 
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iamiltonian 
3, Q, but 



for b —* ssd is obtained in analogy to the one for K° — K° mixing 11, 
with several important differences. First, the CKM structure is more involved. Second, 
the presence of the massive b quark in the initial state introduces a correction, which is 
however suppressed by m^/m^, and is thus numerically negligible. Finally, in applications 
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to K° — K° mixing the charm quark can be integrated out of the theory, while this cannot 
be done for exclusive B decays, where there is no clear separation between the charm mass 
m c and the energy scales relevant in nonleptonic exclusive B decays into two pseudoscalars. 

At scales m b < fi < raw, the effective weak Hamiltonian mediating b — > ssd decays 
contains both local AS = 2 terms as well as nonlocal terms arising from T-products of 
AS = 1 effective weak Hamiltonians 

n ssd = n AS = 2 + jd d xT{nr =1 (x),nt s=1 (o)} . (9) 

The local part is 

n ±s=2 = + xdX(Jct + x , x , Ctc) [{M)v _ A{ , b)v _ A] > (10) 

where the CKM structures are defined as \ q q ' = V qq iV* s . The Wilson coefficient coming from 
the top box loop is C u ~ 0(x t ) and from the top-charm box loop C ct = C tc ~ 0(x c ), where 
Xi = m^/m^y. The scaling of the three contributions in the local Hamiltonian (jTOl) in terms 
of Cabibbo angle A = 0.22 and quark masses is then: ~ X 7 Xt, ~ X 3 x c and ~ X 7 x c (for 
b — > dds all three terms are suppressed by another factor of A). The third term in (fTUl) can 
thus easily be neglected. Note also that there is no A^A* term. The resulting absence of large 
ogx c from the charm box contribution is sometimes called the super-hard GIM mechanism 



14j . and follows from the chiral structure of the weak interaction in the SM, as explained 



in Appendix A. The precise values of the Wilson coefficients in ffTOj) can be read off from 



the expressions for K° — K° mixing llj, where the RG running is performed only down to 
scale (jl ~ m&. For fh t {fn t ) = 160.9 GeV, m c (m c ) = 1.27 GeV, a s (rriz) = 0.118 they are at 
fj, = m b = 4.2 GeV: C tt (m b ) = 1.92, C tc (m b ) = 3.75x c = 9.35 • 10~ 4 at leading order (LO) 
(see appendix [B] for the derivation). 

The nonlocal contributions in b — > ssd transition, Eq. (Q, come from insertions of AS = 1 
effective weak Hamiltonians Ti AS=1 and TL AS=1 (see also Fig. [TJ). The AS" = 1 effective weak 
Hamiltonian TL AS=1 is the same weak Hamiltonian relevant for hadronic B decays 

nt^ 1 = % ( £ ^» E c *<$ - v ^ E C M) . (ii) 

V q,q'=u,c 1=1,2 j=3 

with the tree operators Qf b = (qb) V - A (sq')v-A, Qf >h = (qf3b a ) V -A(s a q'p)v-A, and penguin 
operators ^35 = (sb) v _ A (qq)y TA , Q\ & = (s a b l3 ) v ^ A (q l 3q a ) VTA , where the color indices a,/3 
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are displayed only when the sum is over the fields in different brackets. In the definition of 
the penguin operators Q\_§ a sum over q = {u, d, s, c, b} is implied. The weak Hamiltonian 
Ti^ s=1 follows from (jlip by making the replacement b — > d. 

Using CKM unitarity we can rewrite the CKM factors as V* s V ub = -V t * s V tb - V* s V cb . 
The insertions of tree operators with u and c quarks will generate contributions with CKM 
structure A^A*, that are not present in the local AS = 2 Hamiltonian (flOj) . 

n cc =^f\ d X [ d d x C i C j T{Q^(x)Q^(0)+Q^(x)Q^(0)- 

J i,j=l,2 (12) 

From dimensional analysis, the size of this contribution is roughly 

H cc ~ ^m 2 c X d c \ b c (sd)v- A (sb) V - A , (13) 

which is comparable to (fTOl) and needs to be kept. Another set of contributions of comparable 
size coming from double AS = 1 weak Hamiltonian insertions has CKM structure A^A^. The 
nonlocal contributions proportional to Xf A£, on the other hand, are power suppressed, scaling 
as ml, compared to the corresponding ones in (fill]) , which scale as m\. These contributions 
can be safely neglected. 

The appearance of nonlocal contributions is similar to the situation for K° — K° mixing, 
where the effective Hamiltonian below the charm scale contains the T-product of two AS" = 1 
operators mediating s — > duu transitions, in addition to the local operator {sd)v-A{sd)v-A- 
The only difference is that in exclusive b — > ssd decays the charm quark can not be integrated 
out because of the large momenta of the light mesons in the final state. 

The dominant nonlocal operators have CKM structure AfA^ Eq. (1B8[) and A^A* Eq. (1121) . 
These operators contribute to the physical decay amplitude through rescattering effects 
with DD, Dtt, Dtt, ■ ■ ■ intermediate states. Their matrix elements are suppressed relative 
to those of the top box contribution ~ Cu by X~ A (ml/mf) ~ 30%, which suggests that the 
approximation of neglecting m\ suppressed (but CKM enhanced) nonlocal terms may be a 
reasonable first attempt. 

We leave a complete calculation of the nonlocal contributions for the future and present 
only a partial evaluation of b — > ssd branching ratios by relating the matrix elements of 
the local contributions (flUI) to the already measured charmless two body decays using flavor 
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SU(3). We note that the nonlocal contributions were estimated in Ref. [15( using a hadronic 
saturation model, and were found to be suppressed relative to the local contributions. 

For the purpose of the SU(3) relations to be discussed below, it is useful to rewrite the 
effective Hamiltonian (fTUj) as 

Gp 

Hi = —=V uh V* d KiOi , % = ssd, dds. (14) 
V2 

The operators O t are 

O ssd - = {sb) v ^ A (sd)v- A , O dd s = (db) V -A(ds)v-A, (15) 

and the dimensionless coefficients Ki depend only on the CKM factors and calculable hard 
QCD coefficients. We have 

«W = ^/f W ^ (VuV t lC tt + V cd V; s C ct ), (16) 

V*K) v ub V u d 

and similarly for Kdds- Numerically, the coefficients are (at \x = = 4.2 GeV, with CKM 
elements from {l6| ) 

K ssd - = (6.9 ■ 10- 6 )e i51 ° , ^ds = (1.5 ■ lCT 6 )e-* 74 ° . (17) 

The SU(3) symmetry relations derived below require also the C\ + Ci combination of 
Wilson coefficients, evaluated at the same scale /i = m b . At leading log order this is given 
by 

(C 1 + C 2 )(m b ) = f^MV /23 = r /2 K) = 0.85. (18) 



a s (m h ) 



III. SU(3) PREDICTIONS 



We next show how two body B decay widths for b — > ssd and b — > dds transitions can be 
calculated using flavor SU(3). As a first approximation we neglect the nonlocal charm-top 
contributions, as justified in the previous section. Then the processes b — > ssd and b — > dds 
are mediated in the SM only by the local operators Oi in Eq. (|T5|) . Under flavor SU(3) these 
operators transform as 15 

0„j=15 1/2 , 0^ = 15!, (19) 
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where the subscripts denote the isospin. They belong 
15 in the decomposition of the b — > duu tree operators 

Ci(ub) V -A(du)v-A + C 2 (db) V - A (uu) V -A = 



;o the same SU(3) multiplet as the 
121 



^(Ci + C 2 )(--|i53 /2 



1 IX 



(20) 



These operators contribute to AS" = decays such as B — > 7r7r. The explicit expressions for 
15 operators in (120]) are 



15i/ 2 

153/2 



-L[(bu)(ud) + {bd){uu)\ + iy^[(6 S )(^) + {bd){ss)\ - ^=(bd)(dd), (21) 



— -j=[{bu){ud) + (6tZ)(tZ«)] + —=(bd)(dd) 
v 3 v 3 



(22) 



We list the 6 — > dds, ssd exclusive decays in Table [J for B — > PP and in Table [TT] for 
B — ► PV. The PP final states transform under SU(3) as 1, 8, 27, the operators O ss ^ and 
Odds are in 15, and thus there are only two reduced matrix elements, (8|15|3) , (27|15|3). 
These two reduced matrix elements also appear in the predictions for measured AS" = 
decays mediated by the operators in Eq. (1201) . This means that the B — > PP matrix 
elements of the operators O ss j and Odds can be expressed in terms of AS" = decay 



amplitudes such as A(B C 



7T 1 7T 



and others. A similar analysis applies to B — > PV 



decays, where there are four independent reduced matrix elements of the 15 operators: 
(8s|15|3) , (8^|15|3) , (10|15|3) , (27|15|3). These can again be expressed in terms of phys- 
ical B — ► PV AS = amplitudes. We now derive these relations separately for the B — ► PP 
and B — > PV final states. 



A. B — > PP decays 



We use the formalism of the graphical amplitudes 



which makes the derivation of 



SU(3) decompositions quite intuitive. The two independent reduced matrix elements of the 
15 operator are given in terms of graphical amplitudes l3, Q as 
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10 (27|T5|3> 



1 



C 1 + C 2 



(T + C) 



K 



<t + c), 



ssd 



(8]15|3) 



C 1 + C 2 \5 



-(T + C) + A + E 



K 



' / ^(t + c) + a + e). 



ssd 



(23) 
(24) 
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Transition Mode Amplitude 
b ssd B + -> K+if t + c 

5° -> K°K° t + c 

P s ^/.°tt° 75(« + e) 
B s ^K+tt- -(a + e) 
B s -» K°r) 8 Jl(t + c + a + e) 



dds B + -> K°7T + t + c 

B° ^ K°7r° ^(t + c + a + e) 
B° ^ K-ir + -(a + e) 
P s -» t + c 



TABLE I: B — ► PP exclusive decays mediated by the 6 — > ssd and 6 — > dds transitions. 

This gives two relations between the graphical amplitudes T (tree), C (color-suppressed 
tree), A (annihilation), E (exchange) in the AS = modes (the expression for B — > PP 
decays can be found in [18j) and the corresponding graphical amplitudes t, c, a, e in b — > ssd 
transitions (the decay amplitudes for B — > PP modes in terms of these are collected in 
Tabled]). Equivalent relations apply between AS = and b — *> dds decay amplitudes. 

The most useful for our purposes is the relation fl23|) . This gives the following prediction 
for the exclusive b — > ssd decays 

A(B + -> K + K°) = A(B° -> K°K°) = K * 8 * V2A(B + -> tt+tt ) , (25) 

C\ + C2 

and similarly for the b — > dds decay 

A(B+ -> 7v°tt+) = A(B S -> K°K°) = Kdds V2A(B + -> tt+tt ) . (26) 

Ci + C2 

Neglecting the l/m b suppressed amplitudes e, a one also has 

—A{B S -> iv%) = ^A(5° -> iv°TT°) ~ KflaJ A(£+ -> TT+TT°) . (27) 

2 Kdds Cl + <-^2 

The remaining amplitudes in Table [I] are proportional to e, a. They are l/m& suppressed, 
therefore we do not consider them further. 

The same SU(3) relations hold also for the decays into two vector mesons, B — ► VVKx, 
separately for each helicity amplitude A = 0, ±. For example, the analog of Eq. (1231) is 

A(fi+ - XJ+Jtf ) = A(B° -> ivfKf) = ^L_v^4(5 + - p+p° A ). (28) 
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Transition Mode 



Amplitude 



b -> ssd B+ - 


+ K* + K° 


iy + cy 


B+ - 


+ K + K*° 


ip + Cp 


B° - 


* K*°K° 


ip + £y + cp + cy 


B s - 


-> K*°ir° 


^(ap + ep) 


B s - 


^K°p° 


^(ay + ey) 


B s - 


-> K*° Vs 


yl(tp + cp + ay + ey) 


B s - 




\j\{tv + cy + ap + ep) 


B s - 


» K* + 7T- 


-(op + ep) 


B s - 




-(ay + ey) 


b dds B + - 


-> ^*°7T + 


tp + Cp 


B+ - 




ty + Cy 


B°- 


-» i?*°7T 


73 (*p + c p + + e v) 


B°- 




75 (*V + cy + ap + ep) 


B°- 


» i^*-7r+ 


-(ay + ey) 


B°- 




-(ap + ep) 


B s - 


♦ AT* AT 


ty + tp + cy + cp 



TABLE II: P — > PV exclusive decays mediated by the 6 — > ssd and b — > dds transitions. 

As a consequence the b — > ssd and b — > dds B — > VV decays are longitudinally polarized in 
the same way as the B + — > p + p° decay. 



B. B -> PV decays 



Table HTl lists the decomposition of 5 — > PV" decays in terms of graphical amplitudes. The 
subscripts P, V on t, c identify the final state meson that contains the spectator quark, while 
the subscripts on a, e denote the final state meson containin g th e q 3 quark from b — > q\q 2 q^ 
(here the spectator participates in the weak interaction) 



20 



211 • 



We have T Py + C Py oc (10|15|3) ± (27|15|3). The analogs of the relation ([23]) are then 

^ssd im 1 n \ • , ^ssd 



tp + Cp 



Cx + Q 



(Tp + Cp), t v + c v 



c x + c 2 



(Ty + Cy) 



(29) 
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where the graphical amplitudes on the right-hand side are for AS = decays. The expansion 
of the corresponding decay amplitudes in terms of graphical amplitudes can be found in 
Refs. 
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211 ] . Combining them with expansions in Table HT1 gives the SU(3) relations for the 



ti + Q exclusive b — ► ssd decay amplitudes (for AS* = amplitude we only denote the final 
state) 

A{B + -> K* + K°) = \ - (A p+7T - - A p - n +) - V2A p0w+ 

+ {Ak*°k° — Ak*o k o) — (A K ,- K + — A K »+ K -) + (A^,o^+ — A K *+xo)], 



(30) 



A(B + -> K+iT ) 



(^p+7T- ~ ^p-7r+) — V^Ap+^-O 



c x + c 2 

{Ak*°k° — Ar*°k°) + {Ak*-k+ — Ak*+k-) — {Ak*or+ ~ A-k*+k°) 



A(B U -> ir°K ) 



6(i 



(31) 



(32) 



C X + C 2 

The B s decay amplitudes containing U+Ci are given in terms of the above b — > ssd amplitudes 
A(B S - K*° V8 ) = \j\[A{B + - if+iT ) - A(£? s - KV)], (33) 
- Jjf%) = y|[A(5 + - K*+K°) - A(5 S - iT+O] , (34) 
where the l/mj, suppressed pure annihilation and exchange decay amplitudes are 



A(B S -»■ ir*+7r-) = -v / 2A(fi s -> IT 



7T 



S'6'ci 



^ssd 



[A-k*°k+ ~ A K *- K + — A K *o K o] , (35) 



[A 



K*+K° — A K *+ K 



(36) 



The relations for the b —>■ dds transitions are derived in an analogous way, giving for the 
ti + Q amplitudes 

,*0 A-\ ^dds 



A(B- 



(Ak*°k° ~ Ak*o k o) + (A K *- K + — A K ,+ K -) — (Ak*o k + — A K *+xo) 



A{B + -> i?°p H 



^dds 



C x + C 2 

+ (Aft:* A -0 — Ak»o k o) — [Ajc*-K+ — A K »+ K -) + [Ak*°K+ ~ A k *+ko) 



(37) 



(38) 



and 



V2A(B° 
V2A(B° 



R*%°) = A(B- 
K°p°) = A(B + 



+ iT°7r + ) - A(5° -> K*-7T- 

#V) - a(#° -> . 



(39) 
(40) 



12 



The l/mf, suppressed pure annihilation and exchange amplitudes are 

A(B° -> K*-n+) = [ - A K , +R0 + A K * +K - + A K * oR o] , (41) 

A(B° -> = [ _ ^,o K+ + A K *- K+ + Aj^oxo] • (42) 

The remaining 5 S mode is given by 

A(B S -> = _^L_v^[A p+w o + A p c V+ ]. (43) 

These SU(3) relations will be used in the next Section to predict branching fractions of 
exclusive b —>■ ssd and b — > dds decays in the SM. The measured branching fractions of 
AS" = modes are then the inputs in the predictions and are collected in Table HVl We only 
quote results for those decays that are not l/mj suppressed. 



IV. SM PREDICTIONS FROM THE SU(3) RELATIONS 

Experimentally one will be able to search for NP effects in the following b — > ssd decays 
B° -> K°*K *, B~ K~K *, B° s -> (pK°*. The flavor of K°* is tagged using the decay 
K°(*) _> K + 7i-. The same decays with K° instead of K *, on the other hand, cannot be 
used to probe b — > ssd transitions. The K° mixes with K so that mass eigenstates K$,l 
are observed in the experiment. The "wrong kaon" decays listed above are thus only a 
subleading contribution in the SM rate. For easier comparison with previous calculations in 
the literature we will still quote results for B° — > K°K°, . . . , "branching ratios", knowing 
that these are unobservable in practice. Similar comments apply to b — > dds transitions, 
where NP effects can be probed in B° -> tt° K°* , p° K°* , 5~ -> n-K°*,p-K * and 5° -> 
if * if * decays, again using flavor tagged K°* decays. 

We derive next numerical predictions for the branching fractions of the exclusive b — > 
ssd, dds modes. The branching fraction of a given mode B q — ► M1M2 is given by 

£(5 g -> M X M 2 ) = r B JA(S, -> M!M 2 )| 2 -^L (44) 

To predict b — > ssd, dds decay amplitudes, v4(£> 9 — > MiM 2 ), we use the SU(3) rela- 
tions derived in Sec. II I II which relate them to the amplitudes of the already measured 
B + — > 7i + 7i°,p + p°, and B — > p-n decays. The results are collected in Tables [TTTI and HVl 
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As mentioned, we do not present results for the branching ratios of the l/m& suppressed 
annihilation modes. 

In the calculation of B — > PV branching ratios we neglect the contributions of the small 
penguin dominated B — ► K*K, K*K decays in the SU(3) relations (with experimental upper 
bounds supporting this approximation). Furthermore, the application of the SU(3) relations 
requires that we know also the relative phases of the B — > pn amplitudes. These phases are 
small, and can be neglected to a good approximation. This can be verified using the isospin 
pentagon relation 

A(p + 7T°) + A(p°n + ) = -L(A(p+n-) + A(p-n + )) + V2A(p°n°) . (45) 

V 2 

Neglecting the relative phases, and using data from Table IIV} the left-hand side of this 



equality is 6.25 ± 0.29 (in units of y/B- 10 6 ), which compares well with the right-hand side 
of 6.91 ± 0.34. This justifies the assumption made of neglecting the relative phases of the 
B — > pn amplitudes. 

To factor out the dependence on CKM elements, we also quote the predictions for B 
PP, PV, VV modes in a common form as 

B < B -^(cMr" (46) 

where q are coefficients specific to each final state calculated using the SU(3) relations and 
measured AS* = branching fractions. In the predictions we used the branching fractions 
for the AS = modes listedjn Table El We use t(B + )/t(B°) = 1.071 ± 0.009 and 
t(B° s )/t(B°) = 0.965 ± 0.017 



Both Belle 
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231 ] and BABAR 24( collaborations presented the results of a search for these 



modes and report the 90% C.L. upper bounds (BABAR bounds are in square brackets) 

b^dds: B(B + -> K-it + it + ) < 45.0 [9.5] x 10" 7 ; BELLE [BABAR], (47) 
b^ssd: B{B + -»■ K + K + ir-) < 24.0 [9.5] x 10~ 7 ; BELLE [BABAR], (48) 
B(B° -> K°K + ir-) < 180 x 10" 7 ; BELLE. (49) 

The quasi two-body decay B + — > K°*tt + is part of the B + — > K~7t + tt + three body decay, 
B+ _> K + K*° is part of B+ -> K + K + n~, while B° -> K°K*° is part of B° -»■ K°K + ir~. 
The bounds on three body decays thus imply bound on two-body decays. These are 8 
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Mode 


Ci [xl0- 6 ] 








Literature 


B + 


-> K+K° 


11.0 ± 0.8 


(0.7 ± 0.1) 


10" 


-15 


Z.O X 1U 


-ii 


B° 


-» K -RT° 


10.2 ± 0.7 


(0.7 ± 0.1) 


10" 


-15 






B+ 


— ► 


29.3 ± 4.3 


(1.9 ± 0.3) 


10" 


-15 


1.7 x 10" 


-ii 


B+ 


K+iT* 


11.3 ±3.0 


(0.7 ±0.2) 


io- 


-15 


6.5 x 10- 


-ii 


B° 




71.5 ±6.2 


(4.7 ±0.4) 


io- 


-15 






B+ - 


— ► 


47.2 ±3.7 


(3.1 ±0.2) 


io- 


-15 


6.8 x 10- 


-ii 


B°- 


— > 7^*0^*0 


43.9 ±3.5 


(2.9 ±0.2) 


io- 


-15 







TABLE III: SU(3) predictions for the branching fractions of the b — ► ssd modes in the SM. The 
last column shows the predictions from a previous calculation [3J. 





Mode 


Ci [xl0" 6 ] 


£>SM 




B+ 




11.1 ±0.8 


(35 ± 2) x 10- 


18 


B s 




9.7 ±0.7 


(30 ± 2) x 10- 


18 


B+ 




11.4 ±2.9 


(36 ± 9) x 10- 


18 


B+ 




29.5 ±4.3 


(92 ± 13) x 10" 


-18 


B° 




5.3 ± 1.4 


(17 ±4) x 10- 


18 


B° 




13.7 ±2.0 


(43 ± 6) x 10- 


18 


B s - 


-» K*°K° 


69.1 ±6.0 


(215 ± 19) x 10 


-18 


B+ 


-> tf*V 


47.6 ±3.8 


(148 ± 12) x 10 


-18 


B s - 




41.7 ±3.3 


(130 ± 10) x 10 


-18 



TABLE IV: SU(3) predictions for the branching fractions of the b — ► dds modes in the SM. 

orders of magnitude or more above the estimates for the SM signal, but the situation could 
improve at a future super-B factory 25j or at LHCb. Note that B° — > K°K + tt~ is observed 
in KsK + 7r~ final states which also receives contributions from b — > d penguin decay B° — > 
K°K*° and from annihilation decay B° — K + K~ . It thus cannot be used as a null probe 
of NP. 
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Mode BfxlQ- 



Mode B(xlQ- 6 ) 



B+ - 


-> ir°ir+ 


5.59^ B° - 




23.0 ±2.3 






8.7t\i 


p + 7T~ 


15.4 ± 1.8 a 




p+vr 


W.9±H 


p~7T + 


7.2 ± l.l 6 






2.40 ±0.19 


p°7T 


2.0 ±0.5 



TABLE V: Branching ratios for B — ► 7T7t, p7r, pp decays, from Ref. [22] apart from: a) the average of 



15.5 ±3.4 |26( and 15.3 ±2.2 



, and b) the average of 7.1 ± 1.9 



and 7.3 ± 1.4 



V. b -> ssd AND b (ids TRANSITIONS IN THE PRESENCE OF NP 



Next we consider the 6 — > ssd and 6 — > (ids decays in the presence of generic NP. The 
most general local NP hamiltonian mediating the b —>■ ssd and b — ► dds transitions was 
given in Eq. (TH). In this section we will assume that NP matches onto the local operator 
Qi in Eq. (JTj) with SM chirality (V — A) x (V — A). This is true for a large class of NP 
models, such as the two-Higgs doublet model with small tan /3, or the constrained MSSM 
4j. Effects of NP that matches to other chiral structures will be given in the next section. 1 
To simplify the notation we focus on b — > ssd transitions — the expressions for b — > dds can 
be obtained through a simple s <-> d exchange — but show numerical results for both types 
of decays. 

We consider three representative cases of NP: i) the exchange of NP fields that carry 



in the Introduction, ii) 
NMFV 3 and in ) 



a conserved charge, where large effects are possible as explainer 
minimally flavor violating (MFV) new physics with small tan/5 p 
general flavor violation with a ~ 10 3 TeV scale suppression. For this analysis it is useful to 
rewrite the b —>■ ssd SM effective Hamiltonian (I14D as 

1 



n 



ssd 



G p 

71 



VubV u d K ssdO ssd 



A 2 



e ^ K ssdQl ) 



(50) 



where A = 2 1 / 4 /(2^G F \V ub V ud \) = 2.98 TeV and Qi is defined in © (the flavor dependence 
of Qi is not shown). The NP Hamiltonian for b — > ssd is 

Cl 



n 



NP 
ssd 



A 2 



V2Q1 



(51) 



1 There is also the possibility that NP contributes through the AS = 1 hamiltonians appearing in the 
nonlocal term in Eq. ©, for instance through a (sb)(cc) term. We do not pursue this possibility further. 



16 



where the Wilson coefficient c\ contains possible extra flavor hierarchy in the new physics 
transitions, and Anp is the scale of NP. The hard QCD correction to the c\ coefficient 
describing the RG running from the weak scale m^y to m& has been explicitly factored out, 
V2{ m b) — 0.85. It is now very easy to obtain the branching ratio in the presence of NP from 
the SM predictions, 



B NP (B - /) 
and similarly for b — > dds decays. 



Ao C1V2 



K ssd Anp 



B SM (B - /), 



(52) 



A. NP with conserved charge 

As discussed in the introduction it is possible to have large NP effects in b — > ssd decays, if 
the transition is mediated by NP fields that carry a total conserved charge, and if in addition 
there exists a hierarchy in the couplings. In this case we have for the Wilson coefficient in 
the NP Hamiltonian (ED) (cf. Eq. ®) 

TT~ = ~TFT (db^sQUd + 9d^s9* s -> b ) • (53) 
i\ NP m x 

From K — K and B s — B s mixing we have the bounds, Eq. (pQ), 

\9d^ s gU d \ 1/2 l \9b^ s gUb\ l/2 1 (M) 

M x 10 3 TeV M x 30TeV 1 ' 

These bounds are trivially satisfied, if for instance g s ^d = 9b-^s — 0, since then no mixing 
contributions are induced. The b —>■ ssd transitions, on the other hand, can still be large, 
if gd^s and g s ^b are nonzero. Taking ga^s = gs^b = 1, the BABAR experimental bound 
on B(B+ -> K+K*°), Eq. gHD, gives M x > 5.0 TeV. A similar bound M x > 5.0 TeV is 
found for b — > dds from the BABAR bound on B + — > K*°7r + branching ratio. The resulting 
predictions for b — > ssd and b dds branching ratios are of (9(10 -6 ) as shown in Table [VTl 
and may well be probed at Belle II and LHCb. 

A more generic situation may be that only one of the couplings is accidentally small. 
Unlike in the previous example, we choose M x such that we do not saturate the present 
experimental bounds on b —>■ ssd. As an illustration let us take g s ^d = and all the other 
couplings to be equal to 1. In this case the K — K mixing bound in (15"11) is trivially satisfied, 
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T\ ft 1 / J\ 

Mode (ssd) 


(2 r w /• 5 TeV \4l 

#x[x( ) ] 


r*, / 173 TeV\4l 

lv ssa 


n 




to 

£"SM 




±> — > A A 


1 a i n — 6 
1.4 X 1U 


i n w i n — 12 
l.U X 1U 


n Q w 1 ri- 
ll. O X 1U 


-14 


/'n i _i_ n o\ w i n— 
(U. r ± U.Zj X IU 


•15 


R+ . Z^*+ 1^*0 

— ' i\ J\ 


p. n v i n — 6 

O.U X 1U 


/I 9 v i n — 12 


i /i v 1 n~ 

1 .4t X iu 


-14 


^o.l zt U.Z J X IU 


-15 


D — r J\ I\ 


0.0 X 1U 


q n v i n — 12 
o.y x iu 


i.o x iu 


-14 


l^z.y it u.zj x iu 


-15 


Mode (lids) 


^x[x(^) 4 ] 


^NMPV[X(^) 4 ] 


^gen. 




£>SM 




B+ _» a-*°7t + 


1.4 x 10~ 6 


2.0 x 10~ 14 


1.3 x 10" 


-15 


(3.6 ±0.9) x 10" 


-17 




6.0 x 1(T 6 


8.5 x 10~ 14 


5.4 x 10" 


-15 


(14.8 ± 1.2) x 10" 


-17 


£ s -► A > * A V * 


5.3 x 1(T 6 


7.5 x 10~ 14 


4.7 x 10" 


-15 


(13 ±1) x 10~ 17 



TABLE VI: Predictions for the branching fractions of the b — > ssd, (ids modes in the presence of NP 
carrying conserved charge (2?x)> with NMFV flavor structure (jBnmfv)j an d general flavor violation 
at high scale (£> gen .), in all cases assuming dominance of the SM operator (sd)v-A(sb)v-A (see 
also text for details). The branching fractions B gcn include (maximally constructive) interference 
with the SM amplitude, and were obtained using A ss d t dds = 10 3 TeV. Only modes which do not 
contain Kg^L are shown. 

while B s — B s mixing implies that Mx > 30 TeV. The b — > ssd branching ratios are 
B({B + -> K + K*°,B+^ K* + K*°,B ^ K*°K* }) = {1.1,4.6,4.3} x Kr 9 ( 30l " (A N 



V M 



x 



(55) 

For b — > dds transitions the same choice for the values of coupling, g s ^d = and all the 
other gi = 1, sets Mx > 210 TeV due to the present absence of NP effects in Bd — Bd mixing. 
This gives 



B({B + - K*\ + , B + - K*°p+, B° s - A^A* }) = {0.5, 1.9, 1.7} x 10" 12 ( ^^V- 

V Mx ' 

(56) 

Finally, we mention that, if b — > ssd or b — > dds modes are observed in the near future, 
this would imply nontrivial exclusions on the parameter space of the models. In particular 
models with g s ^ ~ g s ->d and/or g^ s ~ ga-^s would be excluded as discussed in Appendix 



18 



B. NP with MFV and NMFV structures 



Both MFV [28] and NMFV |29J fall in the class of new physics models where the b — > ssd 
suppression scale A ss d is the geometric average of the NP scales in K — K and B s — B s mixing 
(JTJ, A ssd ~ y/ A sd A bs > 173 TeV. In this paper we will restrict ourselves to MFV with small 
tan (3, where the AF = 2 processes are mediated by a single operator with (V — A) x (V — A) 



structure 



28l | . This implies that the K — K and B s — B s mixing operators are 



^^{sliM 2 = tV(^7A) 2 , ^P^{b Ll ,s L f = -^{b Ll ,s L )\ (57) 

MFV 1V sd MFV A fcs 

and the b — > ssd local operator is 

^ 6 y t :^(6 i7 M^)(^7^L) = -t^(&l7^l)(47^l), (58) 



A Z ■ w is- zs\-u ip-uj \-±s l -i-i/ K 

MFV 2 \ s d 



all of which depend only on one unknown parameter, the MFV scale Amfv- From a global 



0. 



fit the UTfit collaboration finds Amfv > 5.5 TeV [9j. We have also defined the suppression 
scales A sd , A bs , A ss j that include the hierarchy of the NP induced flavor changing couplings, 
which in the MFV case are just the appropriate CKM matrix elements. They are related as 



stated above A ssd - = ^A sd A bs . 

In NMFV the operators in ( 1571) and ( 1581) are still parameterically suppressed by the CKM 
matrix elements, but the strict correlation between the Wilson coefficients is lost - they are 



multiplied by 0(1) complex coefficients. We then have approximately A ssd ~ ^A sd A bs . 
Using the bounds from Qj, A ssd > 173 TeV, which gives b — > ssd branching ratios of 
0(1CT 12 ), Table IVTl Similarly we have A dds ~ \/A sd A bd > 458 TeV, giving b — > dds branching 
ratios of O(10~ 13 ). In MFV the predicted branching ratios are much smaller, of the order 
of the SM branching ratios in Table I VII The reason for the difference between NMFV and 
MFV is that in MFV NP the Wilson coefficient generating K — K mixing carries a weak 
phase (the same phase as it does in the SM), while in NMFV the NP contribution can be 
real. 

C. General flavor violation with a high scale 

As a final example we consider the case where NP is at the mass scale probed by K — K 
mixing, A ~ 10 3 TeV, and assume that flavor violating couplings are all of 0(1). The 
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resulting branching fractions for b — > ssd and b — > dds decays, assuming positive interference 
between SM and NP contributions, are collected in Table [VTl For b — > ssd decays the NP and 
SM contributions are roughly of the same size, while for b — > dds the NP induced branching 
ratios are more than two orders of magnitude larger than the SM ones. This means that 
with enough statistics one could probe flavor violation without theoretical uncertainty to 
scales A ~ 10 3 TeV both in 3 — > 2 and 3 — ^ 1 transitions and not just in 2 — * 1 transitions 
as is possible now from K — K mixing. Of course, the statistics needed to achieve such an 
ambitious goal is well beyond the reach of present and planned flavor factories. 



VI. NP LEADING TO NON-SM CHIRALITIES 



We now turn to the description of effects induced by the local operators with non-standard 
chiralities Q2-5 ; Qi-5- It is convenient to normalize the matrix elements of these operators 
to the ones of the SM operator Qi 

and similarly for Q1-5, where the ratio is denoted as fj. To obtain predictions for a b — > ssd 
decay branching ratio due to a particular NP chiral structure, one only needs to multiply 
the results in Table IVT1 with appropriate rj or f|. 

Using parity one can relate r 3 - and fj, since P^QjP = Qj. For B — > PP (B — > VP) 
decays one then has 

fi = =Fl, and fj = =pj > J — 2, ... , 5. (60) 

For B — > VV decays it is convenient to define ratios r\j, f\j for final states with definite 
helicites, (V^aV^a), where A = 0, ±. We then have f lj± = f 10 = —1 and 

fj,± = -r j:T , fj, = -Tj, , j — 2, ... ,5 (61) 

We only need to compute the ratios Tj, j = 2, ... ,5. The ratios f~ are then already 
given by the above relations. To compute Vj we use naive factorization [30j, which suffices 
for the accuracy required here. Strictly speaking, naive factorization is not valid at leading 
order in the heavy quark expansion, but corresponds to assuming dominance of the soft- 
overlap contributions in the complete SCET factorization formula [311 ] . and keeping only 
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terms of leading order in a s (mb). In t 
hard spectator scattering contributions 



re QCDF approach, this corresponds to neglecting 
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341 ] . If needed, these assumptions can be relaxed. 



Naive factorization, or the vacuum insertion approximation, is also justified in the 1/N C 
expansion for the matrix elements of the operators Qi s 2,4, but not for Qs^. To see this, one 
can rewrite Q3 as a sum of color singlet and color octet terms using the color Fierz identity, 

Qs = {s a R bi)(4d a L ) = ^(s R b L )(s R d L ) + 2{s R t a b L )(s R t a d L ), (62) 

1/2 

and analogously for Q 5 . The matrix element of the color-singlet operator scales as N c , 

— 1/2 

while that of the color-octet scales as N c ■ The two term in the above decomposition thus 
contribute at the same order in 1/N C expansion, and both should in principle be kept. 

For the experimentally interesting B — > PV and B —>■ VV decay modes all the ratios can 
be expressed in terms of r 2j 4. One has r 3 5 = 3r 2) 4, and r 4 = — l/[2(iV c + 1)]. 

The ratio r 2 is common to all the PV modes which depend only on the graphical am- 
plitudes tp + cp (for which the spectator quark ends up in the pseudoscalar meson) and is 
given by 

= 1 fyfAmy) g(ro| - m 2 P - m 2 v ) 
V2 8(N c + l)f v f + (m v ) m v (m B + m P ) ' 1 ' 



Using fx* = 218 MeV, /^» = 175 MeV and the form factors from Ref. 35| we find 
r 2 (K + K*°) = 0.28 and r 2 (TT + K*°) = 0.27. 

For the VV modes we quote only the ratios corresponding to longitudinally polarized 
vector mesons, which dominate the total rate. We find rf = —1/[2(N C + 1)] and 

I, _ 1 AmlfT^m^) - {ml ~ m 2 V2 ){ml - m vl - m V2 )T 2 {m vl ) 



8(N C + 1) f vl m V i ( mB + m V2 ){m\ - m Vl - m^AMli) ~ ^&M m " 



.(64) 



m B +m V 2 y v L ' 

Here V\ denotes the neutral K* meson (K*° for b — > ssd transitions, and K*° for the b — > dds 
transitions), if V\, V 2 are different vector mesons. Numerically we find 

r\ (B + -> K* + K*°, K*°K*°) = 0.002 , (65) 
r\{B + -> K*°p + ) = 10~ 5 , r\(B s -> K*°K*°) = 0.002 , 



where we used the B — > V form factors from Ref. 



36]. 
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VII. CONCLUSIONS 



The exclusive rare B decays b — > ssd and b — > dds analyzed in this paper appear in 
the SM only at second order in the weak interactions and have thus very small branching 
fractions, but in NP models they can be greatly enhanced. We construct the complete 
effective Hamiltonian contributing to these modes in the SM, and point out the presence of 
nonlocal contributions, not included in previous work, which can contribute about 30% of 
the local term. 

We show that the hadronic matrix elements of the local operators contributing to these 
exclusive decays in the SM can be determined using SU(3) flavor symmetry in terms of 
measured AS = decay amplitudes. Detailed numerical predictions are given for all B — > 
PP, VP, VV modes of experimental interest, both in the SM and for several examples of 
NP models: NP with conserved global charge, (N)MFV models and general flavor violating 
models. 

A general NP mechanism was identified which can enhance the branching fractions of 
these modes, while obeying existing constraints on NP in AS = 2 mixing processes. This 
mechanism represents a generalization of the sneutrino exchange in R parity violating SUSY. 
Any observation of such a decay mode gives a constraint on the ratio of flavor couplings 
to the NP, and can exclude regions in the parameter space of the NP theory for branching 
fractions observable at LHC-b and super-B factories. 
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APPENDIX A: THE STRUCTURE OF THE EFFECTIVE HAMILTONIAN 

Consider for definiteness the b — ► ssd transitions. At scales M w > fj, > nib, these 
transitions are described by an effective Hamiltonian with propagating u, c quarks. Writing 
explicitly the quarks propagating in the box diagram, and not assuming the unitarity of the 
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CKM matrix, the effective Hamiltonian is given by 

H ssd =\t\ b t H(t,t)+ Aj i A^(g 1 ,g 2 ), (Al) 

with A** = V qq 'V* s (so that for instance \f = V tc iV t * s . The top term in the Hamiltonian is a 
local operator 

n 2 m 2 

H(t,t) = ^-^C(m 2 t /m 2 w ^/m w )[(7sd) v . A (sb) v ^ A ], (A2) 

where C(m 2 /m 2 v , n/mw) is a Wilson coefficient. The box diagram with internal quarks 
Qi,Q2 — u,c, on the other hand, is matched onto an effective Hamiltonian containing both 



local and nonlocal terms 



12| 



m l + ^2)[( sd )v-A(sb) v _ A ] 



H( qi ,q 2 ) = ^ {A(^)m 2 w i(sd) v „ A (sb) v _ A ]+ b(-^)(i 
2 [ \m w J \m w J 

+D^-^)m 2 b [(sd) v ^ A (sb) v . A ]^ + J d d xT{H d (q 2 ,qi)(x),H b (qx,q2)(0)}. (A3) 
The effective Hamiltonian 7~Cb{qi,q2) mediates b — > sqiq2 transitions, and is given by 



ft) = % ( E C *Q$>" + 6 ** E W}) • (A4) 

1=1,2 j=3 

T~id(qi,q2) mediates d scjiq2 and is given by a similar expression, with the replacement 
b — > d. Note that there is no top-charm contribution at leading order in the mf/m^ 
expansion. The top-charm box is matched in the effective theory onto six-quark operators 
of the form (cb) (sc) (sd) , which are power suppressed by relative to the 4-quark 

operators shown. Such terms appear only after using the unitarity of the CKM matrix. 

The dependence on the light quark masses m 12 in the effective theory expression Eq. fl A3j) 
can be obtained in the mass insertion approximation. The W coupling W^(uL"fndL) con- 
serves chirality, which implies that only mf,m| terms are allowed, but not mirri2, which 
would require one mass insertion on each propagating line. The m| term arises from two 
mass insertions on the incoming b quark line. This term is not present in K° — K° mixing. 
On the other hand, in a theory with chiral-odd quark couplings, such as e.g. the charged 
Higgs couplings H + (uLdR) in the 2HDM, another term can appear in Eq. (1A3[) . proportional 
to raim2. Chirality prevents also the appearance of terms of the form mi,mi,m(,m2. 

Under renormalization, the local operator with Wilson coefficient A([i/mw) renormalizes 
multiplicatively, while the nonlocal operators mix into the local operators with coefficients 
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Making use of the unitarity of the CKM matrix, it is possible to eliminate X b u , as 
A^ = —A* — X\, i = b, d. This reproduces the effective Hamiltonian quoted in text Eq. fflOl) . 
The terms proportional to A, C and D are combined into Cu, while the B term in Eq. ( 1A4I) 
reproduces the Ctc and C c t coefficients. The total contribution of the local terms proportional 
to the Wilson coefficient B(fi/mw) is equal to 

■ + K At • ml + A^ ■ ml + X d c X b c ■ 2m 2 c (A5) 
= -X d c (X b c + X b t ) ■ ml - X b (X d + X d ) ■ ml + X d X b c ■ 2m 2 c = -(X d X b + X d c X b )m 2 c . 

This proves the two properties of the local effective Hamiltonian TC AS=2 stated in the text: 
i) the equality Ct c = C c t, and ii) the absence of a X d X^ local term. The latter property 
does not hold in the presence of chiral-odd quark couplings, as for example in the 2HDM as 
discussed above. 



APPENDIX B: AS = 2 WILSON COEFFICIENTS 

In this appendix we show the translation of results obtained for K° — K° mixing to the 



ts for b dds decays are equivalent). The results for 
13 ] in the leading-log approximation, and in [fjjl in the 



case of b — » ssd decays (the resu 
K° — K° mixing were derived in 
next-to-leading log approximation. 

We start with the Wilson coefficient C u , which is obtained by matching the u,c,t loops 
at the weak scale onto the local operator (sb)y_^(sd)y_A- Below this scale, QCD radiative 
corrections introduce a correction ^(/-O, so that at NLO 

C tt (ji) = m(ji)So(xt) + ^^-Dfji/mw). (Bl) 

The box function So(x t ) with x t = m\ /Mir is the same as obtained in the one-loop matching 
at the mw scale for K° — K° mixing (external b quark leg can be considered as massless for 



the purpose of this calculation). It is given by [U] 



Axt-llxl + x 3 , 3x 3 t \ogx t _ 

So[Xt > ~ 4(1 -x t y " 2(1 -x,) 3 " [ } 

with the numerical value given for m t (m t ) = 160.9 GeV. The QCD correction ^(aO is 

obtained by solving the renormalization group equation 

» d = (B3) 
24 



At one-loop order, the anomalous dimension is 7+ = a s /ir, which gives using a>s(mz) = 0.118 
=5 

'MS 
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m(l*b) = [ as ^ Mw } ) = 0.85 , m b = 4.2 GeV, (B4) 



(from which A^T = 226 MeV) 

a, 

a s (m b ) 
so that 

CuifMb) = 1-92. (B5) 

The coefficient D{[x) parameterizes the b quark mass effects, and is introduced by mixing 
from the nonlocal operators into the local operator ml(sb)v-A{sd)v-A- This mixing has not 
been computed yet. We will neglect this contribution since it is suppressed by the small 
ratio ml/rri^ ~ 0.2%. 

The nonlocal contributions due to insertions of two four-quark operators are power 
suppressed and can be neglected as discussed in appendix [A] This is no longer true for 
top-charm contributions, where both local and nonlocal contributions are power suppressed 
by mf/m^r, and mix under renormalization. 

We use the derivation of Jjj], which we adapt to the b — > ssd process at hand. The local 
part of the K° — K° mixing weak Hamiltonian for u above the charm quark mass (i.e. before 

n 

charm quark is integrated out) is given by [14| 

n 2 ~ ~ m 2 
H? S ~ K = ^A?CVQ 7 , Q 7 = -Z[(sd) v -A(8d)v-A). (B6) 

The corresponding local part of the b — > ssd effective Hamiltonian on the other hand is 

G 2 m 2 

-^f(\ d c \ b t C ct + Xt\ b c C tc ) [{5d) V -A(sb)v-A], (B7) 

The RG evolution calculation for b — > ssd process is the same as for K° — K° mixing, except 
that the total contribution is split into two because of two different CKM element structures 
in ( 1B7I) . As shown in Appendix [A], these structures have identical coefficients in the SM 
C c t — Ctc- 

The same equality can be seen also in the anomalous dimension matrices for the running 
of these coefficients. Consider the nonlocal contribution to b — > ssd with insertions of the 
tree operators T{Q 12 Qi^}^ which is given by 

CiCj { x d c x b t {QZdQTb ~ Qiftjfi) + x t x c (Q%Qjfi - QW C A } • ( B8 ) 

25 



When computing the mixing into the local operator Q 7 , the terms in the first and the 
second brackets give the same contributions, since the quark masses are not relevant for the 
calculation of the anomalous dimensions (it does not matter whether c quark or u quark 
runs in the lower leg of the loop in Fig [I]). This shows that the RG running for C c t,Ct c is 
the same. Furthermore, this running is the same as that of C 7 in K° — K° mixing. This can 
be seen by comparing (1B8j) with the nonlocal operator contributing to K° — K° mixing 

£ CiCfAfAj (2QJ3QJ3 - Q%Q% - Q%Qf d ) , (B9 ) 

The two operators are identical, provided that one sets b — > d in ( 1B8I) . The same corre- 
spondence between K° — K° mixing and b —>■ ssd applies also for the nonlocal contributions 
involving penguin operators. 

In conclusion, comparing the Eqs. (IB6I) and (IB7I) we find that for fi > m c , we have 

C ct (fi) = CM = C 7 (fx)x c n/a s , (BIO) 

where C^(fi) is obtained from RG evolution in the same way as for K° — K° mixing. A very 



compact form of RG equations was presented in 

/4-£ = 7 T -A (Bll) 
d\i 

with 7 the 8x8 anomalous dimension, given in Eqs. (6.23)-(6.26) and (12.50)-(12.56) of 



ll| and 

D T = (C T ,C 7+ /C + ,C 7 ./C_). (B12) 

Here C is a vector of Cj, i — 1, . . . 6, C± = C\ ± C2, 2 and C 7 was split to C 7 = C 7+ + CV_, 
where the distribution between C 7+ and C*7_ is arbitrary. At LO we have for the matching at 
weak scale D T (fj,w) = (1, 0, 0, 0, 0, 0, 0, 0), so that the nonzero value of C 7 (/j,) comes entirely 
from the running, from mixing with C\. At \Xb the solution of RG running at LO is 

5 M ^(P^f° V *) V~\ (B13) 

with 7 = 2s. /y(o) anc i y a matrix that diagonalizes the LO anomalous dimension matrix, 
7 g) = rV ^. This gives 

C 7 (m b ) = 0.268, m fe = 4.2 GeV, (B14) 



2 Here we caution about the definition of C1.2, which differs from the one in We use the definition, 
where Ci(/iw) ~ 1, C 2 {^w) ~ 0. 
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and finally 

C tc (m b ) = 3.75x c = 9.35 ■ 1(T 4 , (B15) 
where in the last equality we used m c = 1.27 GeV. 

APPENDIX C: BOUNDS ON THE FLAVOR-CHANGING COUPLINGS 

We have showed in the introduction that b —>■ ssd branching ratios can be large, if NP 
effects are due to exchange of particle(s) with conserved charge. The resulting effective 
weak Hamiltonian, Eq. (j3J, depends on four couplings, g s ^d, 9d^s, 9b^s, 9s^b and an overall 
mass scale Mj, that in this appendix we set to Mx = 10 TeV (this then fixes the overall 
normalization of gi). In order to have large b — > ssd branching ratios and simultaneously 
avoid bounds from K — K mixing and B s — B s mixing a hierarchy between couplings is 
required. Another way of looking at this is that, if a large b — > ssd decay branching ratio 
(we will quantify what "large" means below) is found by Belle II and/or LHCb this would 
imply that a region of parameter space with g s ^ b ~ g s -+d and/or g b ^ s ~ g d ^ s would be 
excluded. We show this below. 

The experimental constraints from K — K mixing and B s — B s mixing give the following 
upper bounds (fixing M x = 10 TeV and using bounds from Eq. ([I])) 

M 2 M 2 

e sd = \g^s9Ud\ < -rr = 10-4 > £ f S = \g^s9Ub\ < -rr = 011 • ( C1 ) 

A sd A bs 

We also define the following two ratios of coupling constants 

R = ^=±^ R=^L. (C2) 

g s ^d gd^s 

We now show that a measured lower bound on the b — ► ssd branching fraction excludes 
values of R, R that are close to 1. For definiteness, we assume that the NP field X couples 
to the quarks with the Dirac structure V = Pr, as in RPV SUSY. Similar bounds can be 
derived for any other Dirac structure T. 

The amplitude for the B — > / transition mediated by the operator (sb)(sd), Eq. (j4j), is 

A {B -»•/) = jp(f\gd-*s9*s^ b Q4 + gb^ s gU d Qi\B) = -^{f\Qx\B){g d ^ s g* s ^ b t gb->s9t+d)> 

(C3) 
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where the upper (lower) sign is for a PP(PV) final state. The combination of couplings gi 
can be written in terms of the ratios R, R defined in (IC2I) 

9b^ s g* s ^ d T gd^gl^b = 9b^ s g* s -+bj£ T 9d^ s g* s ^ d R* = {g d ^ s g* s -> d )R T (<fo-»»<£U)^ • (C4) 

The products of coefficients on the r.h.s are now exactly the ones bounded from the meson 
mixing, Eq. (IClj) . The absolute value of the l.h.s on the other hand is assumed to be 
bounded from below from the measurement of b —>■ ssd branching ratio, cf. Eq. (IC3I) . We 
then have 

B 2 < \g b ^ s gU d =F g d ^ s gU b \ 2 < s 2 sd \R\ 2 + 4 77^ + ^s d e bs . (C5) 
If B > 2yje sd e bs , then the above inequality rules out a range of values for \R\, 

^{B- 2e sd e bs - ^/B 2 - Ae sd e bs ] < \R\ 2 < —[B - 2e sd e bs + B 2 - 4e sd e bs ] . (C6) 

Z£ sd Z£ sd 

The same bound with e sd <-> e bs holds also for |.R|. The requirement B > 2^/e sd e bs corre- 
sponds to the requirement that B(B — > /) > 4B(B — > /) NMFV , with the NMFV predictions 
for branching ratios given in Table IVI1 
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